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Abstract: Mutually unbiased bases generalize the X , Y and Z qubit bases. They 
possess numerous applications in quantum information science. It is well-known that 
in prime power dimensions N = (with p prime and m a positive integer) there 
exists a maximal set ofN+1 mutually unbiased bases. In the present paper, we derive 
an explicit expression for those bases, in terms of the (operations of the) associated 
finite field (Galois division ring) of N elements. This expression is shown to be 
equivalent to the expressions previously obtained by Ivanovic in odd prime dimensions 
(J. Phys. A, 14, 3241 (1981)), and Wootters and Fields (Ann. Phys. 191, 363 
(1989)) in odd prime power dimensions. In even prime power dimensions, we derive 
a new explicit expression for the mutually unbiased bases. The new ingredients of 
our approach are, basically, the following: we provide a simple expression of the 
generalised Pauli group in terms of the additive characters of the field and we derive 
an exact groupal composition law inside the elements of the commuting subsets of 
the generalised Pauli group, renormalised by a well-chosen phase-factor. 



Introduction 



Two orthonormal bases of a dimensional Hilbert space are said to be mutually 
unbiased if whenever we choose one state in the first basis, and a second state in 
the second basis, the modulus squared of their in-product is equal to It is 
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well-known that, when the dimension of the Hilbert space is a prime power, there 
exists a set of + 1 mutually unbiased bases. This set is maximal because it is not 
possible to find more than + 1 mutually unbiased bases in a iV dimensional Hilbert 
space [HElinj. It is also a complete set because when we know all the probabihties of 
transition of a given quantum state towards the states of the bases of this set (there 
are N"^ — 1 of them), we can reconstruct all the coefficients of the density matrix that 
characterizes this state; in other words we can perform full tomography or complete 
quantum state determination PUSS El- A crucial element of the construction is the 
existence of a finite commutative division ring (or field^) of elements. As it is well 
known, finite fields with N elements exist if and only if the dimension is a power 
of a prime, and a derivation of a set of mutually unbiased bases is already known in 
such cases. Note that nobody managed until now to generalize this construction in 
the absence of finite field so that it is still an open question whether such sets exist 
when the dimension is not a prime power [HIIZI. In the present paper, we obtain, 
in a synthetic formulation, the expressions for the mutually unbiased bases that 
were derived in the past. In odd prime power dimensions, we recover by a slightly 
different approach the expressions already obtained in the past by Ivanovic and in 
odd prime power dimension by Wootters and Fields (2j). We provide a synthetic 
expression that is also valid in even prime powers dimensions (2"^). The (discrete) 
Heisenberg-Weyl group jSl IH E] (sometimes also called generalized Pauli group), a 
finite group of unitary transformations, plays a central role in our approach. 



1 Preliminary concepts 

In what follows, we shall systematically assume that we work in a Hilbert space of 
prime power dimension = with p a prime number, and m a positive integer. 
Then, as is well known, it is possible to find a finite field with A^ elements. We 
shall label these elements by an integer number z,0<i<A^ — l),or, equivalently, 
by a muple of integer numbers (zq, ii, im-i) running from to p — 1 that we 
get from the p-ary expansion of i: i = Y^^=q inP^- This field is characterized by 
two operations, a multiplication and an addition, that we shall denote Qg and (Bg 
respectively. It is always possible to label the elements of the field in such a way that 
the addition is equivalent with the addition modulo p componentwise. As all the 

field is a set with a multiplication and an addition operation which satisfy the usual rules, 
associativity and commutativity of both operations, the distributive law, existence of an additive 
identity and a multiplicative identity 1, additive inverses, and multiplicative inverses for every 
element, excepted. 
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fields are equivalent, up to a relabelling, there is no strict obligation to do so, but it is 
more natural and convenient. This property is a direct consequence of the fact that 
for all the finite fields the characteristics of the field, which is the smallest number 
of times that we must add the element 1 (neutral for the multiplication)with itself 
before we obtain (neutral for the addition), is always equal to a prime number {p 
when = p"*). The index G refers to Evariste Galois and is introduced in order not 
to confuse these operations with the usual (complex) multiplication and addition for 
which no index is written. 

Let us denote 7g the pth root of unity: = e''^'^/^. Exponentiating jg with 
elements g of the field (with the usual rules for exponientation), we obtain complex 
phasors of the type 'Jq {0 < g < N). Such phasors can take p different values. They 
can be considered as ap-uple generalisation of the (binary) parity operation e*'^^'^/^^'^ 
that corresponds to the qubit case in the sense that the phasor 'Jg {0 < g < N) only 
depends on the value of the first component go of the p-ary expansion of g which is 
nothing else than the remainder of g after division by p, when the division by p is 
taken in the usual sense. 

The following identity appears to play a fundamental role in our approach: 

N-l 

E 7r"^^ = A^^.o (1) 

j=0 

Indeed, if ^ = 0, then Ef=o' 7^®'''^ = N.1 = N. Otherwise, Efjo tJ^'®'''^ = E?=o 7g 
in virtue of the inversibility of the multiplication. Now the exponentiation of gamma 
by elements of the field does only depend on the remainder after division by p, so 

max 2^ji=Q 1g —V ■ 2^j'=o 1g — V ■ (1-70) — ^■ 

In virtue of the fact that the addition is the addition modulo p, componentwise, 
we can derive the following identity which is also very useful: 



lh-lh = 7^^^'^ (2) 

indeed, Tg ' Tg ~ ig ih = Xg = 7g l^^ previous expression, 
we represented by the symbol xq the remainder of x after division by p, where x is 
an element of the field, comprised between and p"* — 1 = — 1, and the division 
by p is taken in the usual sense.) This relation is well-known and expresses, in the 
language of mathematicians, that pth roots of unity are additive characters of the 
Galois field 
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It is important to note, in order to avoid confusions, that different types of 
operations are present at this leveh the internal field operations are labelled by 
the lower index G. They must not be confused with the modulo operations. In 
order to emphasise the difference between these operations, we give in example the 
corresponding tables in the case iV = 4 = 2^ in appendix. One can check that 
the field and modulo 4 multiplications are distributive relatively to the associated 
addition, but that there are no dividers of 0, excepted, only in the case of the 
field multiplication. As a consequence, the field multiplication table, amputed from 
the first line and column exhibits an invertible (group) structure. All operations 
are commutative as can be seen from the symmetry of the tables 1 to 4 under 
transposition. 

Remark that if we express quartits as products of two qubits: |0)4 = |0)2 |0)2, 
|1)4 = 10)2 ® 11)2, 12)4 = 11)2 ® 10)2, 13)4 = 11)2 ® |1)2- It is then easy to check the 
following property: If |ii)4 = \ii)2 ® 1^2)2, and \j)4 = \ji)2 ® ^2)2, then \i ©g j)4 = 
1^1 ©mod2 Ji)2 ® 1^2 ©mod2 ^2)2- TMs is an illustratiou of the fact that the field addition 
is equivalent with the addition modulo p componentwise. It is also worth noting 
that the property J2p=o 7^^®'^-' = NSq^ is true for the modulo N multiplication as 
well, but 7 must be taken to be equal to the A^th root of unity in this case. In prime 
dimensions 7^ is the A^th root of unity and the Galois and modulo operations 
coincide. In prime power but non-prime dimensions, this is no longer true. 



2 Construction of the dual basis 



Let us now consider the unitary transformations VJ^, that shift each label of the 
states of the computational basis ({|0), |A^ — 1)}) by a distance / — >• 

\i ©G 0) (the reason for our choice of notation will be made obvious soon). The 
transformations Vf form a commutative group with A^ elements that is isomorphic 
to the Galois addition. Generalizing the procedure outlined in 0, we define the 
dual basis as follows: 

1 N-l 

\k) (3) 

where the symbol ©g represents the inverse of the Galois addition ®g- It is easy to 
check that the dual states are invariant, up to a global phase, under the transfor- 
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mations VJ^. Indeed, we have: 



1 N-l 

^Mj) = 4^E7S^^'^^^'^I^©gO (4) 

■1 N-l 

= 4^E7r^''^^'^^^^'^|fc')=7r^^'^b) (5) 

Obviously, the dual basis and the computational basis are mutually unbiased. When 
the dimension is prime {N = p), the dual basis is the discrete Fourier transform of 
the computational basis, when it is a power of 2, it is a Hadamard transform p. 

Let us denote Vq the unitary transformations that shift each label of the states of 
the dual basis ({|0), — 1)}) by a distance QqI — >■ \i Qg 0)- The 

transformations Vq form a commutative group with elements that is isomorphic 
to the Galois addition. It is easy to check that these operators are diagonal in the 
computational basis: 

7V-1 _ _ _ N-l 

K=T.\kQGl)Ck\=T.^G'"''^\k){k\. (6) 

A:=0 A,-0 

This is the dual counterpart of a similar expression for the shifts in the computational 
basis: 

N-l N-l 

^°=EI^©G0(fc| = E7r"'^l^)(^l- (7) 

A;=0 k=0 

3 Construction of the remaining mutually 
unbiased bases 

In the previous section we derived a set of two bases, the computational basis and 
the dual basis, that are mutually unbiased. In this section, we shall generahze 
this derivation in order to obtain other mutually unbiased bases (between each 
other, and also relatively to the computational and dual bases). 

Let us denote the compositions of the shifts in the computational and the 
dual basis: 
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N-l 

= VIV:^ = E 4''''''^'''''^\k®G^){k\■,^,J : 0...N - 1 (8) 

fc=0 

Here, the product . expresses the matricial product (the usual composition law of 
two unitary transformations). 

As Vq is the identity, there is no disagreement with the previous definitions. 
Note that Vq and Vj' do not commute: 

N-l N-l 

V/ = vivr = E V®^^IO(/|- E \k ®G ^){k\ 

1=0 k=0 
N-l 

fc=0 

Although this expression appeared, to our knowledge, in ref.[19^ for the first time, 
we show in the last section that the set of operators so-defined coincides with the 
generalised Pauli group studied in ref.jH]. 



N-l N-l 

k=0 1=0 
N-l 

k=0 

^ ^eG(«0Gi)T/i T/.° = -v©G(«0Gi)T/; 

/ ■ ' i I ^ i 

The commutator is thus given by the following expression: 
Vi.V^ - V^.Vi = (1 - ^^oi^<^o^)^vlVl 

We recognize here a commutation rule that is known as the Weyl commutation 
rule, and was already studied a long time ago 0. This is not astonishing because 
the set of unitary transformations V!' that we consider here is a discrete version of 
the so-called Heisenberg-Weyl group (compositions of translations in position and in 
impulsion). In dimension 2, it coincides with the Pauli group. When the dimension 
is a prime number, the field operations are the addition and multiplication modulo p, 
and the properties of mutually unbiased bases are already well-known in that case , 
as well as their relation with the "Heisenberg-Weyl-Pauli" group lOj. In the present 
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(10) 
(11) 



approach, we consider, instead of the usual (modulo A^) operations, the Galois 
addition and multiplication, also for non-prime but prime power dimensions. The 
connection with previous works related to the Pauli group approach [3 El HB 1121 CHI 
is established in the last section. 

By a straightforward computation, we can now derive the law of composition of 
these iV^ unitary transformations: 

= (12) 

Up to a global phase, this looks like a groupal composition law. We shall now 
show that (up to phases) the N"^ unitary transformations form iV + 1 commuting 
subgroups of elements that have only the identity in common. Moreover, each of 
these subgroups admits a diagonal representation in a basis that is mutually unbiased 
relatively to the the N bases in which the other subgroups are diagonal. Note that 
the last property can be shown, following an alternative approach developed in ref. 
12] to be a consequence of the fact that the V operators, up to phases, form what 
is called a maximally commuting basis of orthogonal unitary matrices (see also last 
section). The new ingredients in our approach are (1) the expression El for the 
generalized Pauli operators, and (2) the recognition of the fact that these operators 
exhibit exact groupal composition laws, provided (a) they commute and (2) they 
are multiplied by a convenient phase factor. 

In order to derive all the results, we shall take them for granted in a first time, 
and check afterwards that our hypothesis was correct. It is convenient to introduce 
new notations and definitions before we pursue. We shall denote Ul the elements 
of these subgroups, where i labels the subgroup and runs from to (there are 
+ 1 of them), while I labels the elements of the subgroup and runs from to A^ — 1 
(each subgroup contains A^ elements). We know already the two first subgroups, that 
admit a diagonal representation in the computational and dual bases: the first one 
{i = 0) contains the elements V^{1 : O...A^-l), so by definition Uj^ = 1/J(/ : O...A^-l). 
The second one contains the elements Ul- = VP{1 : O...A^ — 1). 

In virtue of the equalities El and we can also write f/° = S£o^ 7g^'''^I^)(^I 
and ul = J2k=o 1g'^'^^'^\^){M- ^ similar expression can be found for each of the 
A^ — 1 remaining subgroups as we shall now show. It is convenient at this level to 
parametrize the basis states that diagonalize these subgroups as follows: the fcth 
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basis state that diagonalizes the ith subgroup will be denoted as |e^). Our ultimate 
goal is to prove that there exist 1 bases |e^) and N"^ operators Ul that are in one 
to one correspondence with the V operators and differ from them by an appropriate 
phase factor, such that the following constraints are fulfilled: 

t^/= E7r"'^l4)(e^l (/:0...iV-l;.:0...iV) (13) 

k=0 

(ef |e;.).(e5|ef) = 6kA,j + {^/N).{1 - {kj: O...N,t,j : 0...N - 1) (14) 

In virtue of the commutativity of the Galois multiplication, of the identity El and 
of the definition El the U operations that are labelled by a same value i form a 
commutative subgroup and obey the (exact) group composition law Ul-Up = Uli^^i,. 
We can guess that they correspond to families of operators V/' such that the (Galois) 
ratio k/c^ is constant, because the commutation of V^'^ and V^^ implies that k'Qcl = 
k Qg I' ■ It is thus natural to try the identification f// = v^^*"-'^)®^;'^ ^ phase, 

when i differs from and = Vq which is consistent with our previous conventions. 
There are in general several ways to fix the phases but in any case certain constraints 
must be satisfied: 

-the phase f///V^'^*"^'*®'^' is equal to when / = 0, because Vq = 1, and the 
identity is present in all subgroups 

-as we mentioned already, the U operators must obey the following composition 
law: Ui^^i, = Upl,, but the composition law F/.V^*^ = ^QG{i®Gk)Yi®Gk ^^^^ 
guaranteed at the same time, which restricts seriously the arbitrariness in the choice 
of the phase. 

Let us now assume that the phase ratio between Ul and Vi'' is fixed for all 
powers of p between and m — 1 {I = p", < < m — 1). 

We shall firstly treat the odd dimensional case. Then, iterating I times the 
composition law^l (with 2 < Z < m — 1), we obtain the following constraints on 
the ratio between U^n^^i and V^t~^^f''^"®''\ 0</<m-l,0<n<m-l: 

N-1 
k=0 

— /TP /T/(*-1)QgP"V ^eG(i-l)0G«0G(«eGl)0GP"0GP"/G2 t A({i-1)0GP"0GO ( ^ ,-\ 

— [Up„/ Vpn ) ■ iG ' >"0G^ V-*-"^/ 



8 



Here, the symbols / and / g indicate the multiplication by the multiplicative inverse 
for the usual (complex) and field multiplications respectively. In order to fix the 

phase (t/^^®^^"V^p^'^®''''") 

we can make use of the fact that the characteristics 
of the field is p (so to say 1 ©g 1 ®g ■■■ ®g 1 (p times) = 0), which implies that 
(C/pn)^ = 1, so that we obtain the following constraint: 

The phase {Ul/Vpn ^^^'^-^ ) is thus determined up to an integer power of 7g = e'"^'^/^. 
This is true for each integer value of n between and m — 1 so that there are — N 
different possible ways to "fix" the phases. Let us denote 7„ the pth root of unity 
that we choose to be equal to Upn/Vpl . Once this value is chosen, all the 

other phases are determined, as shows the following development: 

tH _ T-rm-lTTi _ rrm- 1 / r ri \ln 

— ■'■^n=0 <^;„OgP" ~ ■'-^n=0 V<^p"J 

_ TTm-l^,eG(j-l)0Gin0G(«neGl)OGP"OGP"/G2/ , V„ T/((i-l)0GP"0Gin) /lyN 
— ■lin=0 7G Kin) >"0Gi„ ) K^l) 

where the coefficients /„ are unambiguously defined by the p-ary expansion of /: 
I = Ek=o InP"'- Moreover, we can check by direct computation that the U operators 
so-defined obey an exact groupal composition law, independently on the choice that 
we could decide to perform, among the — N different possible ways to "fix" the 
phases 7„: 



Tji Tji _ TT'"~^7'r TP 

— TT"»-lo,©G(i-l)0G'ln0G(^ln©Gl)0GP"0GP"/G2„0G(i-l)0G'2„0G('2„©Gl)0GP"0GP"-/G2 

— ■ii„=o 7g 7g 

\hn(^ V2nT/((i-l)OGP"0G;in)T/((j-l)0GP"0G«2n) 
V In) V '^p"0G'ln '^P"0g'2„ 

— TT"l-l„,©G(i-l)0GP"0GP"0G('ln0G('lneGl)©G'2n0G('2neGl)©G20G'ln0G'2„)/G2 

— iin=0 7g 

(r^ Vl„+modp'2nT/((«-l)QGP"0G(iln+moV2n)) 

V 'n) ^ P^Qoihu+modphn) 

_ TTm-l„ ©G{«-l)0Gj'"0GP"0G{('lTi+modp'2n)0G({'ln+modp'2n)©Gl)/G2 

— iin=0 7g 

\h„+rr,odpl2n-\/((^-'^)®GP"QGihn+modpl2n)) 

V ^ P''QG{hn+modphn) 

- iin=0 - ^^n=0 l^(ii„+Wp'2„)0GP"^ ~ ^/i©Gi2 \^^) 

In even prime power dimensions, the treatment is similar, although we may not 
divide by 2 in this case. Combining the constraints (^1!^"^^®°^"^ = (^eGO-i)0G2"0G2" 
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7g and UlnUln = t^2"eG2" ~ = 1; obtain the following decompo- 

sition law for the U operators, which expresses their factorisation in terms of qubit 
operators: 

'^1 — -'-'-n=0 '^«„0g2" ~ -'-'■n=0 2^ ) 

where the coefficients In are unambiguously defined by the binary expansion of I: 
I = jy^=o ^n2"; In = or In = 1- Similary to what happens in the odd dimensional 
case, the phase factors (7^"^''^'^^ ^'^'^ )^ can be fixed with some arbitrariness, 
actually up to a sign in this case. Let us now check that the U operators so defined 
obey an exact group composition law. 

U^-.^^t'^'^-""''''"")-'''^^ = UUair (20) 

We made use of the fact that = ^O-i)0g2"0g2"^ ^ 

Let us now derive an explicit expression for the iVphases Ul/Vi^~^'^®°\ We shall 
treat separately even and odd prime power dimensions. 



3.1 Odd prime power dimensions 



In odd prime power dimensions, all possible consistent choices for determining the 
phases (there are = N such choices) can be expressed as follows: 

f^i/^(i-l)0G« ^ ^^e((i-l)0G«0GO/G2^^^fc0i^ ^21) 

where k is an arbitrary element of the field. Each choice for k (there are of them) 
leads to another consistent determination of the phase ratio between the U and V 
operators. This is due to the fact that if jq^'^p'^ = Vn : < n < m — 1, then 

k = k' in virtue of the identity H 

Note that, when k = 0, which is the simplest determination of the phases, we 
obtain the following relation: 

Ui = (^g©((»-l)®G;0GO/G2^y^(»-l)0Gi ^22) 
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This corresponds to the choice of phases 7„ = ry'^^'y''-~^)®GP @gp )/g2^ j^. .^qj-^j-^ 
noting that the relation |221 is also valid for i = 1, which corresponds to the dual 
basis derived in the previous section. 

For i = 0, Uj^ = Vq, in agreement with the previous definitions. 

In order to check the consistency of the expression fEgnl^Tj). it is sufficient, 
making use of the composition law for the V operators fEqn J12|) . to check by direct 
computation that the U operators obey an exact (so to say not up to a phase) group 
composition law. Remark that if Ul = 70^^®'?//, and that f//.f/// = U^^^i,, then 
U}.U}, = U}^^i,. Therefore, it is sufficient to establish the groupal composition law 
when the expression fEqnl22j) is valid, in order to establish it when the expression 
fEqnJ21|) is valid, for any value of k. 

Ui.Ui, = ('y^©((»-l)QGteGO/G2^(^^^e((i-l)0Gi'0G«')/G2^^(*-l)0G«^(i-l)0G«' 
= (^^p(^-'^)®GlQGl)/G2^^^^e({i-l)QGl'eGn/G'2^^^(eG(i-l)QGl^ 

= (7g®«*-')®«('®«'')®«('®^''^)/«')V^(S;"^\^^^ = fZ/ecF;^ : l...N,l,l' : 0...iV- 1(23) 

Now that we have at our disposal an exact expression for the operators U, we 
can also derive an explicit expression for the — 1 dual bases associated to the 
subgroups that correspond to the operators Ul;i — 1 = l...A^ — 1. This can be 
realised thanks to the following identity, a direct consequence of EqnsHni and [H : 

1 N-l 

l4)(eII = ^E7®^'=®^^t// (24) 
1=0 

Obviously, if we choose another determination of the phases, so to say if we replace 
Ui by Ui = 'yc'^'^^Ui, we obtain the same basis states, with their labels shifted by k'. 
It is thus more convenient to choose in what follows the simplest phase determination 

M 

By a straightforward but lengthy computation that we do not reproduce here, 
we obtain then the expression, in the computational basis, of the states of — 1 
bases that correspond to the non-null values of the label i — 1. 

1 N-l 

14) = ^ E 7|^''®^'(7G«^-^)®-''®-'^)/-^)|e;), (25) 
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Remark that the previous expression is also vahd when i = 1, which corresponds 
to the dual basis \j) 1^. 

Let us now check by direct computation that the bases (A^ — 1 plus the dual 
basis) obtained so are orthonormal and mutually unbiased between each other (it is 
easy to check that the computational basis also fulfills these requirements). Before 
we do so, we shall rewrite the factors (^(^Q((*~^)®G90Gg)/G2-) follows: 

(^^^e({i-i)0G«0GO/G2^ = f/i/V^^*"^^®''' = (^©((*-i)®G'0GO^f ^2Q) 

This redefinition is less precise than the previous one, because there exist two deter- 
minations of the square root of a complex number, nevertheless we adopt it, having 
in mind that in odd prime power dimensions, the previous expression fixes the sign 
of the square root of (^®((*~^)®<3^0gO^ without ambiguity. We shall show in the next 
section that a similar expression is also valid in the even dimensional case. Let us 
now prove that the expression El is valid. 

(e'.|e^) = — 7|''^®'^^*®'^-^'^(7^^''~"^^®'^^'""^^^®'^^®'^''')^ (27) 

^ 9=0 

1 N-l N-1 

(Yl 7^'^''®°^^®°-'''' (7^^^""^^®''^^""^^^®'^'^®'^'^)^) (Y 7®°'^'®°^-''®''*^(7^'~-^^®'^*^''"-^^^®'"''®°''')i) 

5=0 q'=0 

1 Af-1 7V-1 
^ _(^ ^0Gq0G{»©Gi)^^«fc-l)©G{^-l))0G'?0G9^i^ ^0G (q©Gt)0G {i©G«) 1)0G (fc- 1))0G {5©Gt)0G faeG*) ^ i ^ 

q=0 t=0 

1 N-l 

^ — -(^ ^*0G(i©Gi)(^^^(a-l)0G(fc-l))0Gg0Gt^|i j^^a«-l)eG(fc-l))0Gt0Gt^)i-j 

q,t=0 
1 Af-1 

_( ST ^,((*QG(i©Gi)©G({'-l)©G('=-l))0G90Gt) ^^,(('- 1)©G (fc- l))0Gt0Gt> 



^Ut0GU©G«J©GU«-ij©Gi's-iJj0G90Gtj ^^(^ 

q,t=0 



N-l N-l 



1 ( ^'( ^,a('-l)©G(fc-l))0Gi0Gg)-^ ^(t0G(i©Gi)) (^^a'-l)©G (fc-l))0Gt0Gt^ i ' 



t=0 q=0 



1 N-l 

^ sr , , , jt&Gijeai)) / (a-i)©G(fe-i))0Gt0Gta 



1 

r A-LHA ^ 'S:^A o,(*®GO'©Gi)) /^,((«-l)©G(fe-l))0Gi0GtNi 

= 5fcA. + (W)-(i-M(2 
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We made use of the fact that there is no divider of excepted itself (the mul- 
tiphcation Qg forms a division ring). Henceforth, the following identity is valid: 

Finally, let us control the validity of the postulated expression IT^ 



N-l 
k=0 

N-l 1 N-l N-l 

fc=0 ^ q=0 q'=0 

AT ^-1 

q,q'=0 

N-l 

- sr /Ai-'~'^)®G(q'®Gi)@G(q'®Gi)eii-i)QGq'QGq')\h^o \/^o i 

q'=0 

N-l 

= (7c«^-^)««'®-'))^ E 7;,"^-'^®"'^®"^'^|eJeG^)(e;,| = : 0...N - l;z : l...iV)(29) 

(j'=0 



3.2 Even prime power dimensions 

In this case the explicit expressions for the mutually unbiased bases are less easy to 
manipulate. Once again, there are (2™ in this case) possible ways to determine 
the phases U{/Vi ^ i)0g'^ -^^^ they are equivalent, up to a relabelling of the basis 
states. 

In the next development, we shall implicitly choose a certain determination of 
the square root of 7^"^^®^^"®''^" that is equal to i(i"i)0G2"0G2"_ 

— ■'-■'-n=0 ^i„0G2" ~ ■'-'■n=0 l^2"J 

= n;ri.,„(7g-"^'°"^"=")i(i^i'-'"""')'" 

- n'"-! „-{i-l)0G2"0G2'VT/((i-l)0G2")N 

where the coefficients In are unambiguously defined by the p-ary (here binary) ex- 
pansion of /: / = jy^=o ^n2", while n' is the smallest integer strictly larger than n 
such that In' 7^ 0, if it exists, otherwise. 
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This result is a generalisation of the identity !^ because in both cases the phases 
are square roots of integer powers of gamma: 

Nevertheless, in the present case we obtain the following determination of the 
square root of (7g©(0-i)QGteGO): 

(^^^e{(i-i)0G'0GO)5 = (7^©({i-i)0GteGO)| 

= n--i„,o^(-^)««-«--7r^^^^"^^"' (31) 

where n' and were defined previously. What is particular with even prime powers 
is that the square root of an integer power of 7g is not an integer power of ■jg as 
was the case in odd prime power dimensions. We are forced to introduce +i and -i. 
What is interesting is that this minimal extension is sufficient in order to diagonalize 
the operators of the generalized Pauli group in even prime power dimensions, a fact 
that was already recognized in previous references on the subject f^l ITH]). 

Combining the equations ESI and EHl we obtain the following synthetic expression: 



1 N-l 

4) = ^7 E 7r''®^'^(7G«^-^)®«''®«''¥|e;) (32) 



N ,=0 

Taking account of the Eqnj^we get an explicit expression for the mutually unbiased 
bases: 

1 N-l , 

\el) = ^7 E TS^^^^'^n-- „^o^(-^)«-^"«-^"7r^^^''"®^" (33) 

viV q=o 

where q = Y^^=q Qn'2"', while n' is the smallest integer strictly larger than n such 
that Qn' 7^ 0, if it exists, otherwise. 

As a consequence of the composition laws [T21 and EDI 

TtJ Tjj ^ ( U-i)QGheGh\^ / O-l)0G'20G«2a O-l)0G{«i0Gi2)T/(i-l)0GaieG'2) 

h h \ iG J ■{ iG ) I 'i©G^2 

_ /^/i-l)0Gai©G'2)0Gai©G«2)a T/(i-l)0G('l©G'2) _ Tjj 

- \1g I -^h^Gh -'^h®Gh y'^^) 

Therefore, 

^^O-l)0G'lWG'l^^_^^U--Lj0G«2WG«2^^_^(j-lj0G(«10G'2) ^ (^^i 



/j-l)0G^10G'l^|i ^^,(i-l)0G^20G^2^|i ^^{j-l)0G{h0G/2) ^ ^^O■-l)0Gal©G/2)0G(/l©G^2)^| 1 



(35) 
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Formally we can rewrite the previous equation as follows: (7^ i)0G(a©G'')0G(affiGfe)^ ^ _ 
^^a-i)0G{a0Ga))i_(^O-i)0G{;>0Gb))i_(^2.(O-i)0Ga0Gfc))i^ ^j^i^j^ is reminisceut of the equa- 
tionEl although we are dealing here with half integer powers of jg instead of integer 
powers. Thanks to this property, it is possible to reproduce nearly litterally the 
proofs given in odd prime power dimensions of the validity of the identities ^1 and 
El because the automatisms of computation are nearly equivalent. It is impor- 
tant to note however that in even prime power dimensions the expressions of the 
type (7g^'^"'')^ do well represent square roots of 7^®'^"-' but must be considered 
as functions that depend on the 2"^ variables a instead of only two variables, as 
would be the case if we considered litterally square roots of integer powers of 7^ 
(with 7^ = —1). When a is specified, the sign of the square root is also specified, 
according to the explicit expression |^ The even and odd dimensional cases are 
covered by the synthetic expression |221 



4 Open questions, comments and conclusions. 
4.1 Other symmetries 

At first sight, the computational basis plays a special role in our approach, but one 
can show that, to some extent, all the mutually unbiased bases can be treated on 
the same footing. This can be seen as follows. Now that we have at our disposal an 
explicit expression fEqns 03'l33p for all the mutually unbiased bases, we can "reeval- 
uate the situation from the point of view of one of them" , say the ith basis (with 
i different from zero). In order to do so, we can express the action of the operator 
in terms of its basis states. After a straightforward computation, we get that 
V-{0) = phase.V-^^^^,_,^^^Jz), where K^(0) = ^£0^2'®""^^"^"^^ |e^eG..)(e°l 

and V;'J(i) = E£oSG''®''"'^®'"'^l4eGm)(4l5^ • T hese relations (that we give 

without proof but are easy to derive from Eqns ]25in^ are bijective. So the whole 
discrete Heisenberg-Weyl group is invariant (up to permutations and phase shifts) 
when we reexpress it in any of the N +1 mutually unbiased bases. We shall not de- 
velop this question here, but this property has important implications in the theory 
of cloning machines, in relation with error operators and optimal cloning ([SIE])- In 
prime dimensions, the invariance of the Heisenberg-Weyl group under conjugation 
by any unitary matrix that maps the computational basis onto a mutually unbiased 
basis is a basic property of a larger group that is known as the Jacobi or Clifford 
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group and possesses many applications in number theory and quantum computing 

m 



Beside, there exists a one to one correspondence between generahsed Bell states 
[TH] and the Heisenberg-Weyl group. The properties of invariance of the Bell states 
in mutually unbiased bases appeared to be very useful in the resolution of the so- 
called mean king problem jTUl UHl UII , where it also led to a compact and elegant 
expression valid in all prime power dimensions |18j . 

4.2 Connection with previous works. 

The expression of the states of the mutually unbiased bases that we derived in the 
present paper (|ei) = ^ EJTo' 7S'''®°*'(7g is actually equivalent 

to the solution derived by Ivanovic P when the dimension is an odd prime (which 
can be shown, when rewritten according to our conventions, to be equivalent to the 
expression |el)^™"- = ^ Ef=o' tS'^'^'^'^ITg "'^®'"'®'"'^)|e°)). Our expression differs 
by a factor 1/g2 in one of the exponents of 70- When the dimension is prime and 
odd, it is easy to compensate the difference by a relabelling of the basis states, 
because the division by 2 is a permutation of the finite fields with p elements when p 
is a prime odd number, but contrary to Ivanovic's expression, our expression is easily 
generalized in even prime dimension 2 (the qubit case), in which case we rederive 
the eigen bases of the Pauli operators, and in prime power dimensions. 

As it was shown by Wootters and Fields j2], the generalisation in prime power 
dimensions of Ivanovic's expression is the following: 

1 N-l 

\el) = ^ E 7S^-^^-«^'=^(7S^-^^"^^"^^^)|e;), (36) 

where Tr. represents the field theoretical trace. Our expression seems to be a bit 
simpler, but both expressions require to know the addition and multiplication tables 
of the field, so that the apparent gain in simplicity is relative. Moreover, both 
expressions are equivalent up to a relabelling in odd prime power dimensions as we 
shall now show. We could establish the equivalence at once but we prefer to base 
our derivation on the results of the reference [Hj where the interrelation between 
the Pauli group approach and the expression of Wootters and Fields with the trace 
factor is made (section 4.3., jS]). Actually, there exist also different groups that 
present properties similar to those of the generalized Pauli group [T^ but do not 
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obey the definition that we gave here. Nevertheless, the generators of the two 
subgroups corresponding to the computational and the dual basis that are given in 
the reference coincide with our choice (the subgroups Vq and correspond to the 
classes Cq and Ci studied in pj) so that we are talking exactly about the same group, 
up to phases. As, in the same paper, a relation was established with the solution of 
Wootters and Fields [2], our expression for mutually unbiased bases must necessarily 
coincide with the expression derived by Wootters and Fields. In the reference |.3], it 
is shown that when there exists a maximal commuting basis of orthogonal unitary 
matrices, the + 1 bases that diagonalize these classes are unabiguously defined 
and, moreover, are mutually unbiased. A maximal commuting basis of orthogonal 
unitary matrices is a set of A^ + 1 sets of A^ — 1 commuting unitary operators (or 
classes) plus the identity such that these A^^ operators are orthogonal regarding the 
in-product induced by the (usual operator) trace denoted tr.. It is easy to show that 
the V operators defined in Eqn|Hlare unitary with (V/)'^ = (V^/)~^ = ^sgC^Ogj)^^^* 
and that tr.V/ = N.6i^Q.6jfi. Making use of the composition law^l we obtain the 
relation tr.((V^"')'*".V^'^) = N.6i^i.6j^k so that they form a maximal commuting basis of 
unitary operators. This theorem suggests another way to derive an expression for 
the mutually unbiased bases: it is sufficient to find the common eigenstates of the 
classes of operators VJ*'* ^''®'^' (where / varies from to A^ — 1) in order to determine 
the value of the states of the ith mutually unbiased basis. When the dimension is 
an odd (even) prime power, one can check by direct substitution of the expression 
EH (jSSl) that the states |e^) are common eigenstates of the ith class: 



Af-l 1 N-l 

k'=0 V ^ g=0 

^ N-l 



^«'?®GO0G(i-l)0GO^©Gg0Gfc(^^^((i-l)0Gg0Gg)/G2^|gO^ 



q=0 
1 N-l 

^ ^g0Gfc)^({i-l)0G«0GO/G2)_^ J2 ^eG('3®GO0Gfc(^^^((i-l)0G(g©GO0G(g©GO)/G2^|gO^ ^^^37^ 

Thanks to the product law ESI the proof is entirely similar in even prime power 
dimensions. 

In order to establish explicitly and once for all the equivalence between the 
expression of Wootters and Fields (j36p and ours (|25p. some work remains to be 
done because in the reference jS] no proof is given of the fact that the expression 
represents eigenstates of the generalised Pauli operators. In order to prove this 
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result, it is useful to introduce two (field theoretical) dual bases: the first one, 
which is dual relatively to the trace contains m elements p of the field such that 
Tr.p* QgP' = the second one which is dual relatively to the rest after division 
by p contains m elements of the field such that (p* ©c p-^^ = where (g)o 
represents (we work in dimension p™) the rest after division of q by p. These bases 
can be shown to exist and to be unique, in virtue of the fact that the bilinear 
forms Tr.(x ©g y) and (x ©g y)o are non-singular 02!, a direct consequence of the 
identity H Beside, the m generators of the fcth class considered in [3] are equal 
to X*.n]^o^(n^^(Z-^')''»j with j,l : 0...m - 1,. In the previous expression, the 
coefficients ki are unambiguously defined by the p-ary expansion of fc: k = I]/=o^ ^iP^ 
while the multiplication matrix h is defined as follows: 7^ ©g 7j = Ylb=o^ ^^i^^' where 
the 7's are a basis of the field (here we shall consider without loss of generality 
that 7j = Pi). The operators are "local" operators that shift the ith component 
of the label of the kth. basis state |e°) (with k = Ya^^ hp^) by unity (modulo p): 
X*|e°) = |e°/) with k[ = ki +mod.p 1, ^ = ki when I i. The operator multiplies 
the state |e°) by a global phase equal to 7J. In our approach, the generators of the 
k'th class can be shown to be the same, provided the coefficients of k' expressed in 
the double-tilded dual basis p^ defined here above are the same as those of k in the 
direct basis p: k = I]]=o^ ^iP'' ^' — ^T=o^ kjpK Beside, the expressions with and 
without Trace I2H1 and inni are equivalent, up to a bijective relabelling, in virtue of the 
following identity: Tr.(r ©g k) = {{r' /g'^) ©g ^)o with k and r arbitrary elements of 
the Galois field with N = p^ elements, r = E/=o^ '^W^ r' = ^'iP^ ©g 2. 

This comparison emphasises the difference between our approach and previous 
approaches: our expression |H1 of the generalised Pauli operators is global and non- 
local, although they can be decomposed as products of local operators. It also shows 
that the field theoretical trace is replaced in our approach by another non-singular 
bilinear form, the rest after division by p. 

Although it is out of the scope of the present paper, it would be interesting to 
understand the relation between our results in even prime power dimensions and 
the results presented in references [21 US]- 

4.3 Other dimensions. 

It is still an open question to know whether maximal sets of mutually unbiased 
bases exist in arbitrary dimensions. For instance in dimension 6 which is the small- 
est dimension that is not a power of a prime, nobody knows whether or not such 
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a maximal set exists [HI El- It is not possible to apply our treatment in this case 
because no finite field with 6 elements exists. We could try to repeat the procedure 
with operations that do not form a field; for instance we could try to find a dis- 
tributive ring with 6 elements (such a ring obeys the same definition as a field (the 
definition that was given at the beginning of the paper), excepted that the multipli- 
cation needs not be invertible-dividers of zero different from zero are allowed). One 
can show that there is only one distributive ring with 6 elements, that corresponds to 
the usual operations (multiplication and addition modulo 6). If we study the struc- 
ture of the N"^ = 36 Heisenberg-Weyl unitary transformations in that case, we find 
that there are more than + 1 = 7 subgroups of 6 elements (5+the identity). This 
is because, as a consequence of the non-invertibility of the multiplication modulo 
6 (3 and 2 divide zero), certain operators present degeneracies and belong simulta- 
neously to different subgroups (a treatment of similar type is given in detail in the 
reference ^H] for the case = 4 ). The bases that diagonalize these operators are 
not mutually unbiased in general and the construction that was succesfuUy apphed 
in prime power dimensions does not provide a maximal set of mutually unbiased 
bases. Therefore the question of the existence of 7 mutually unbiased bases in a 
6 dimensional Hilbert space is still open, and our approach does unfortunately not 
contribute to the elucidation of that problem. 

4.4 Conclusions 

As we already mentioned, there is a one to one correspondance between (gener- 
alised) Bell states and (generalised) Pauli operators [121 (^^^ ^^^^ ^| for a different 
approach based on additive and multiplicative characters of the Galois field). It can 
also be proven ITU] that the Bell states are invariant when we pass from one of 
the mutually unbiased bases to another one, an important result in the theory of 
cloning machines that was only conjectured until now jHj. Actually, the present re- 
sults were largely inspired by results that we obtained in the framework of quantum 
cryptography P where the interest of mutually unbiased has been recognised several 
years ago, for what concerns encryption [2011211122] and cloning as well [3 [111123 • 

It is worth noting that, beside quantum cryptography and quantum cloning, 
the Bell states found also many applications in quantum teleportation and dense 
coding and the connections between mutually unbiased bases, complete orthogonal 
families of unitary matrices, and teleportation, were already emphasised in the past 
|241 [23] . There exists also an impressive litterature about the interrelation between 
finite fields and discrete Wigner representation [23 [23 [23 • It is worth noting that 
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if we perform a tomographic development of the density matrix in the basis of the 
V operators, we obtain [19, the following identity: p = (l/N) Y.k,i Vi^Tr-^iyi^Y .p). 

It is very instructive to compare the amplitudes of the decomposition of a density 
matrix into the V basis with the Wigner function: 

The Wigner function can be written as follows in terms of the conjugate 
continuous variables q and p\ 

W{(1-,P) = C. J d^rp{q — r, qj^r)e^-''-^-''/^ = C. J d^r'p{r', 2.g — r')e~*-^-P'''/'' where 
C is a normalisation constant, while i is the square root of -1. The analogy of both 
expressions is striking and, as we can see, the tomography of a quantum state that we 
realize in the Pauli group approach provides a discrete counterpart of the Wigner 
representation. In general the coefficients Tr.((VJ^)"'".p) are complex, which does 
not meet the requirements of a properly discretized Wigner function, but in even 
prime power dimensions this difficulty can be overcomen because the V operators 
are Hermitian up to a global phase that was defined in Eqn|^ (the U operators 
defined by EqnJHOIare Hermitian and unitary and also provide an orthogonal basis). 
It is out of the scope of the present paper but it would be interesting to study the 
connection with other proposals for discrete phase-space representation |2Z1 123 12^] • 
Note that as the V operators are diagonal in the + 1 mutually unbiased bases, 
full tomography can be obtained by performing + 1 von Neumann measurements, 
as was already shown by Ivanovic in prime dimensions and Wootters and Fields in 
prime power dimensions. 

Finally, the properties of Bell states are also directly related to the error opera- 
tors [HOI 1^ Hn] , and it would be worth investigating to which extent our formalism 
contributes to a simplification of the theory of error correcting codes, in prime power 
dimensions. 

To conclude, we note that, despite of the fact that the problem (and its solutions) 
seem to be regularly rediscovered by different generations of physicists, which means 
also a lack of time and energy, our results about the Mean King's problem [TH] 
confirm that it is important to explore alternative approaches in the treatment of 
the question of mutually unbiased bases. 

We wish that the present paper contributes to a deeper understanding of the 
old problem of mutually unbiased bases. 
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Table 1: The field multiplication in dimension 4. 
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Table 2: The field addition in dimension 4. 
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Table 3: The multiplication modulo 4. 
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Table 4: The addition modulo 4. 
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